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ABSTRACT 

In  order  to  investigate  the  effects  of  gear  train  backlash  on  the 
performance  of  second  order  servomechanisms,  an  automatic  digital  com- 
puter was  utilized  to  obtain  solutions  of  numerous  second  order  servo 
problems.  Solutions  were  obtained  for  variations  of  motor  and  load 
inertias,  motor  and  load  frictions,  combined  system  damping  coefficient, 
and  backlash  angle,  throughout  the  practical  range  of  these  parameters. 
The  data  generated  by  the  digital  computer  solution  has  been  analyzed 
and  presented  in  graphical  form  for  u^e  in  estimating  the  perforraance 
of  an  arbitrary  second  order  servo  system  with  backlash. 

The  program  was  conducted  on  the  Control  Data  Corporation  Model 
1604  Computer  located  at  the  U.  S,  Naval  Postgraduate  School,  Monterey, 
California, 

The  writers  wish  to  express  their  appreciation  for  the  assistance 
and  encourageiuent  given  them  by  Professor  George  J.  Thaler  of  the  U,  S. 
Naval  Postgraduate  School » 
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Chapter  I 

INTRODUCTION 

A  servomechanism  which  is  designed  using  linear  theory  rarely 
perforiES  as  estimated.  The  discrepancies  between  design  performance 
and  actual  performance  are  caused,  in  general,  by  certain  non-line- 
arities which  are  inherent  in  all  systeiTiS.  Of  the  many  non-line- 
arities which  are  found  in  servo  systems,  the  one  under  investigation 
here  is  backlash. 

In  a  servo  system,  whenever  torques  are  transmitted  from  a  motor 
to  a  load  through  a  gear  train,  there  will  exist  a  certain  amount  of 
backlash  betv;een  the  gears  as  shown  on  figure  1,  Present  day  gear 
drives  cannot  transmit  loads  continuously  without  some  backlash.  If 
the  gears  are  too  closely  meshed  they  might  bind,  wear  excessively,  or 
introduce  random  torque  variations  appearing  as  noise.  This  is  partic- 
ularly true  cf   gear  trains  which  operate  over  wide  temperature  varia- 
tions. Backlash  is  often  controlled  so  as  to  introduce  &   slight  dither 
to  a  system  to  provide  good  lubrication  to  overcome  stictional  effects. 

It  is  generally  known  that  backlash  in  a  servo  system  is 
destabilizing.  However,  the  estimation  of  the  effect  of  backlash  on 
the  transient  response  and  the  ultimate  steady  state  operation  of  a 
system  presents  a  very  difficult  problem.  As  soon  as  backlash  is 
introduced  into  an  otherwise  linear  system,  response  is  then  dependent 
upon  the  division  of  inertia  between  motor  and  load,  the  division  of 
viscous  friction  between  motor  and  load,  the  damping  coefficient  of  the 
system  when  operating  as  a  unit,  the  amount  of  energy  dissipated  as 
heat  when  contact  is  made  between  the  gear  teeth,  and  other  phenomena 
which  ea:-.not  be  predicit-d.  Since  the  system  performance  is  affected 
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differently  in  the  presence  of  backlash  when  a  different  set  of  variables 
is  used J  each  system  presents  a  new  problem. 

In  practice^,  knowledge  of  the  performance  of  a  system  may  be  desir= 
able  prior  to  the  construction  of  a  mock-up  or  simulation  system.  Ji^ 
order  to  obviate  the  necessity  of  estimating  the  performance  by  labori- 
ous and  complicated  analytical  techniques,  design  charts  are  presented 
in  this  thesis  from  which  the  design  engineer  can  predict  the  magnitude 
of  maximum  overshoot  for  any  arbitrary  second  order  servo  with  backlash. 
A  qualitative  analysis  of  the  effect  of  backlash  on  the  tijne  of  maximum 
overshoot  and  the  settling  time  is  also  presented,^ 

Previous  work  has  been  done  in  this  area  at  the  Naval  Postgraduate 
School  by  W.  J.  Lutkenhouse,  (Ref.  l),  and  N.  C.  New,  (Ref.  2). 
Lutkenhouse  solved  the  problem  by  construction  of  the  phase  trajectory 
on  the  phase  plane.  New  extended  the  work  of  Lutkenhouse  by  employing 
an  automatic  digital  computer  for  the  solution  of  the  problem,  and  was 
able  to  collect  a  large  voluirie  of  data  for  use  in  the  construction  of 
design  charts  for  predicting  the  size  of  limit  cycles  in  the  system. 
The  work  done  by  Lutkenhouse  and  New  was  limited  to  the  existence 
and  size  of  limit  cycles  in  a  system.  This  thesis  extends  that  work 
to  include  transient  behavior  of  load  motion.  The  computer  program 
written  by  New  was  utilized  for  the  solution  of  a  large  number  of 
problems^  the  data  from  which  was  subsequently  analyzed  and  processed 
for  the  construction  of  design  curves. 


Chapter  II 
BACKGROUND 

Prior  to  the  development  of  high  speed  digital  computers,  it  v/as 
necessary  for  the  engineer  to  use  long  and  laborious  methods  for 
analyzing  the  behavior  of  servomechanism  systems  containing  non-linear 
elements.  Previous  investigation  of  methods  used  to  analyze  non-linear 
systems  has  resulted  in  several  approaches. 

One  of  the  most  vddely  used  methods  for  analyzing  the  behavior  of 
a  non-linearity  is  the  phase  plane  diagram.  At  present,  this  method 
is  limited  to  control  systems  represented  by  a  second  order  differential 
equation  in  which  the  coefficients  can  be  constants  or  functions  of  the 
dependent  variable.  The  solution  to  this  problem  is  one  in  which  time 
is  not  expressed  explicitly,  but  only  the  velocity  auid  position  of  the 
dependent  variable  for  each  instant  of  time  are  showr^.  From  Eef .  1, 
it  can  be  noted  how  tedious,  time  consigning,  and  often  inaccurate  such 
a  method  is  for  obtaining  the  transient  response  of  a  second  order 
servo  with  some  non-linearity. 

Another  method  widely  used  in  non-linear  problems  is  the  use  of 
describing  functions.  This  method  uses  the  concept  of  obtaining  an 
"equivalent"  transfer  function.  If  a  sinusoidal  input  of  constant 
amplitude  is  applied  to  a  non-linear  system,  the  output  will  contain 
a  number  of  harmonics  of  the  input  frequency.  If  the  fundamental 
component  of  the  output  is  assuj-ned  to  be  of  the  same  frequency  as  the 
input,  then  an  "equivalent"  transfer  function  can  be  obtained  which 
is  the  ratio  of  the  fundamental  component  of  the  output  to  the  sinusoi- 
dal wave  input.  Once  this  relationship  of  the  output  tc  the  input 
of  a  non-linear  system  is  obtained,  it  is  then  used  in  the  same 


manner  as  the  transfer  function  in  the  linear  systemo  Since  this 
method  assumes  that  the  higher  harmonics  have  negligible  effect  on  the 
system  and  that  the  frequency  of  the  input  is  the  same  as  the  output, 
it  is  seen  that  it  can  only  be  used  for  a  rough  approximation  in 
determining  the  amplitude  and  period  of  the  output.  This  method  is 
also  disadvantageous  because  of  the  difficulty  in  the  computations 
associated  with  the  determination  of  the  describing  function  for  a 
wide  range  of  amplitudes  and  frequencies.  There  is  no  simple  method 
for  evaluating  the  accuracy  of  a  describing  function  and  no  good  way 
to  check  the  validity  of  solutions  obtained  with  it. 

With  the  advent  of  the  analogue  computer,  another  method  of 
calculating  system  transient  performance  v;as  available.  This  method 
merely  simulates  the  problem  of  the  servo  and  its  npn-linearities  to 
obtain  the  response  of  the  system.  This  method  is  of  good  value  in 
the  preliminary  evaluation  of  the  system  behavior  where  accuracy  is 
not  too  important.  However,  due  to  the  inherent  inaccuracies  of  the 
analogue  computer  this  method  does  not  give  precise  results. 

By  utilizing  the  digital  computer  in  the  solution  of  non-'linear 
servos^ it  is  possible  to  overcome  many  of  the  disadvantages  inherent 
in  the  other  methods.  A  non-linear  system  which  can  ba  described  by 
a  few  accurate  or  approximate  differential  equations  can  be  solved 
on  a  digital  computer  resulting  in  transient  and  steady  state  informa- 
tion. In  this  thesis  a  digital  computer  was  used  to  solve  many  second 
order  servo  problems  with  backlash  present.  From  the  wealth  of  results 
obtained,  the  response  of  any  arbitrary  second  order  servo  with  back- 
lash may  be  predicted.  It  is  believed  that  this  procedure  could  be 
applied  to  servos  of  higher  order  and  with  other  non-linearities. 


Chapter  III 
DISCUSSION  OF  SYSTEM  BEHAVIOR 

The  system  under  investigation  in  this  thesis  is  a  simple  second 
order  position  servo  control  having  the  motor  coupled  to  the  load  through 
gear  trains  having  backlash.  The  configuration  of  the  system  is  illus- 
trated in  block  diagram  form  in  figure  2,  Position  feedback  is  taken 
from  the  load  shaft  including  the  backlash  vdthin  the  feedback  loop. 
When  backlash  is  outside  the  feedback  loop,  i.e.,  position  feedback 
from  the  motor  shaft,  the  system  is  stable  and  does  not  result  in  a 
limit  cyclei  however,  it  is  subject  to  a  residual  steady  state  error 
in  response  to  a  step  input  which  error  may  be  as  great  as  the  magnitude 
of  the  backlash.  When  included  in  the  feedback  loop,  backlash  is  de- 
stabilizing and  the  effect  is  dependent  upon  many  variables.  When  a 
step  position  input  is  imposed  on  the  system  there  is  initially  a  period 
when  the  backlash  is  taken  up,  and  the  system  performs  as  a  linear  sys= 
tern.  If  the  linear  system  damping  coefficient  is  less  than  the  critical 
value,  the  error  signal  will  eventually  cause  the  motor  shaft  to  reverse 
its  direction  of  rotation,  A  separation  between  the  motor  and  the  load 
will  occur  in  the  cycle  of  motion  whenever  the  deceleration  of  the  motor 
becomes  greater  than  that  of  the  load.  However,  if  the  frictional 
forces  on  the  load  produce  a  deceleration  greater  than  the  maximum  de- 
celeration of  the  motor,  the  motor  and  the  load  will  maintain  -  ontact 
during  this  period  and  become  separated  only  when  the  motor  velocity 
reverses.  The  precise  instant  at  which  motor  and  load  separate  depends 
upon  the  nature  of  the  frictional  forces  acting  and  the  division  of 
inertia  between  the  motor  and  load,  l^en  the  gears  do  break  contact, 
the  load  drifts  and  the  system  is  no  longer  a  closed  loop  system  but  is 
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driven  open  loop.  Whenever  the  total  amount  of  backlash  is  taken  up 
and  gear  contact  is  made,  the  action  which  follo\vs  is  very   complex. 
The  basic  laws  which  govern  the  ensuing  motion  are  the  law  of  conserva- 
tion of  energy  and  the  law  of  conservation  of  momentum.  In  this  thesis 
the  assumption  is  made  that  the  contact  is  a  plastic  one  with  no  gear 
bounce,  and  that  energy  is  dissipated  as  heat  due  to  the  plastic  contact. 
Thus  only  the  law  of  conservation  of  momentum  are  applied.   Following 
the  gear  contact  phase,  the  systein  is  again  combined,  and  the  cycle  of 
events  begins  again. 

The  motion  of  the  motor  and  load  is  illustrated  in  figiires  3  and  4. 
Angular  velocity  versus  angular  position  is  plotted  and  the  phase  tra- 
jectories for  both  motor  and  load  are  shown  for  two  conditions  of  inertia 
and  friction  ratios.  In  figure  3  it  is  noted  that  at  point  A,  the  sys- 
tem reacts  to  a  step  input  of  1  radiitn  by  moving  in  a  trajectory  in 
accordance  with  the  differential  equation  for  the  linear  system.  At 
point  B  the  gears  which  couple  the  motor  to  the  load  separate 5  the 
dotted  trajectory  from  B  to  C  showing  the  path  followed  by  the  motor 
and  the  solid  line  from  B  to  C  showing  the  path  of  the  drifting  load. 
The  gears  separate  whenever  the  velocity  of  the  load  exceeds  that  of 
the  motor,  which  depends  upon  the  inertia  and  friction  distribution 
between  motor  and  load.  In  this  case  the  inertia  of  the  motor  is  equal 
to  the  inertia  of  the  load,  but  the  load  friction  is  much  greater  than 
the  motor  friction.  For  this  reason  the  gears  do  not  separate  until  the 
motor  has  sufficient  error  signal  in  the  opposite  direction  to  cause  the 

A  study  is  being  carried  on  concurrently  to  include  the  conserva- 
tion  of  energy  and  the  elastic  impact  case  by  Lt.  N,  0,  Anderson  and 
Lt.  T.  W.  Luckett  at  the  U,  S.  Naval  Postgraduate  School, 
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motor  to  decelerate  faster  than  the  load.  As  the  error  signal   (v/hich 
is  determined  by  the  load  position)  increases,  the  ruotor  velocity 
decreases  to  zero  and  begins  accelerating  in  the  opposite  direction, 
V/hen  the  load  and  motor  positions  differ  by  an  amount  equal  to  the  back- 
lash angle  as  sho^wn  at  point  G  and  G*,  the  gears  make  contact  and  the 
lav/  of  momentum  is  applied  causing  the  motor  and  load  to  immediately 
assiime  the  same  velocity^  that  velocity  being  dependent  upon  the  rela- 
tive masses,  or  inertias j,  of  the  motor  and  load.  This  new  velocity  is 
shov^  at  points  D  and  D  ,  The  motor  and  load  now  both  move  in  accord- 
ance v.dth  the  combined  m.otion  equation  but  in  the  opposite  direction 
and  separated  by  an  amoxmt  equal  to  the  Loagnitude  of  the  backlash,  as 
shown  from  D'  to  S*  and  from  D  to  E. 

The  same  cycle  of  events  occurs  throughout  the  motion  of  the  system. 
The  effect  on  the  first  overslioot  is  observed  from  this  much  of  the  tra- 
jectory. 

Figure  4  shov/s  the  ef  f e  rt  of  backlash  on  a  system  v/hich  has  a  load 
inertia  greater  than  the  motor  inertia  ind  a  motor  friction  greater  than 
the  load  friction,  with  the  same  backlash  as  in  the  first  example.  It 
is  important  to  notice  at  point  B  that  the  gears  separate  earlier  than 
in  the  previous  example.  Both  the  heavier  load  inertia  and  lighter  load 
friction  contribute  to  its  slower  deceleration.  As  before,  v.hen  the 
motor  and  load  are  separated  by  the  amount  of  the  backlash  the  gears 
mai<:e  contact  and  the  m.otor  and  load  assunie  the  same  velocity  as  shovm 
by  point  D  and  D",  Mote  that  the  combined  velocity  is  positive  in 
this  case  because  of  the  high  load  moment-uni.  Thus  the  combined  motion 
is  in  the  positive  direction  with  the  error  signal  causing  the  motor  to 
act  as  a  brake  on  the  system  motion. 
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These  two  cases  serve  to  indicate  the  fact  that  every  case  nust  be 
studleci  individually  since  every  coiabination  of  inertias  and  frictions 
cause  considei'ably  different  results. 

In  order  to  make  the  problem  mathematically  tractible,  the  follow- 
ing asi.ijnptions  have  been  made; 

1,  Backlash  is  the  only  non-linearity  present  in  the  system. 

2,  Backlash  is  present  only  in  the  gears  coupling  the  motor  to 
the  load, 

3,  The  gear  contact  is  plastic  thus  applying  the  lav;  of  conserva-= 
tion  of  momentuni  only. 

It  is  believed  that  the  results  obtained  with  these  assumptions 
give  a  good  approximation  of  the  effects  of  backlash  on  a  second  order 
servo  system. 


Chapter  IV 
SOLUTION  OF  EQUATIONS  OF  MOTION 

ViHien  operating  as  a  combined  system  the  servo  motion  can  be 
described  by  the  differential  equation. 

When  the  load  drifts  free  within  the  backlash  region  the  load  motion 
follows  the  differential  equation, 

and  the  equation  for  the  motion  of  the  motor  is. 

The  characteristics  of  the  second  order  servo  with  backlash  is 
independent  of  the  imdamped  natural  frequency.  The  proof  is  develop= 
ed  in  Ref.  2.  Therefore,  all  of  the  problems  solved  and  used  for  the 
construction  of  the  design  charts  are  based  on  an  undamped  natural  fre- 
quency of  one  radian  per  second.  The  input  for  all  problems  was  a  step 
input  of  one  radian.  Thus  equation  (1)  becomes, 

The  differential  equations  IV-=1  ,  17-2,  and  IV-3>  for  the  second 
order  servo  with  backlash  were  solved  on  the  Control  Data  Corporation 
Model  1604  Computer,  The  program  used  was  based  on  the  Rimge-Kutta  Gill 
numerical  integration  method  of  solving  differential  equations.  A  com- 
plete description  of  the  program  and  an  explanation  of  the  solution  are 
given  in  Ref,  2.  The  program  is  included  in  Appendix  C.  By  merely  insert- 
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ing  into  specific  memory  locations  the  system  parameters  for  the  problem 

to  be  solvedj  the  solution  could  be  obtained  for  any  desired  second  order 

servo, 

A  study  of  the  describing  differential  equations  reveals  that  there 

are  five  system  parameters  which  are  independent.  These  are  4y*  ^m'  *^'* 

Fj.,  and  A  .  To  facilitate  the  graphing  of  the  results,  ratios  of 

J       F 
L  and   L    were  selected,  and  from  these  selected  ratios  discrete 

''m     ^total 

values  of  the  parameters  J^ ,  F, ,  J,,,  and  F,,,  were  calculated  for  use  in 

the  con^uter  solution. 

The  system  damping  coefficient  is  determined  by  the  relationship. 


^r  = 


ror^/> 


Taking  J_,  =  J   +  Jj^  =  1  simplified  calculations,  but  did  not 
destroy  the  generality.  The  system  problem  was  solved  for  all  combina- 
tions of  the  following  ranges  of  inertia  and  friction  ratios  and  damping 
coefficients 

^  =  1/50,  1/9,  1/U,   1,  4,  9,  and  90 

M 

Ft 
j~—  -  0,1,  0,2,  0,4,  0,6,  0,8,  1.0 

total 

J^       =  0.1,  0.2,  0.3,  0,4,  0,5,  0,6,  0,8,  1.0 

The  manner  in  which  the  results  are  presented  in  the  final  design 
charts  pennits  interpolation  through  the  practical  range  of  system 
parameters. 

The  amount  of  backlash,  A  ,  used  in  the  majority  of  the  problems 
was  0,3  radian  and  it  was  a  purpose  of  this  thesis  to  determine  the 
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effect  of  other  magnitudes  of  backlash.  At  the  outset  several  problems 
were  solved  varying  the  amount  of  backlash  while  holding  the  other  para- 
meters at  constant  values.  Figures  5-7  are  plots  of  the  results  of  these 
solutions.  Two  items  of  significant  interest  were  noted  from  these  re- 
sults. First,  the  amplitude  of  maximum  overshoot  varies  linearly  with 
the  magnitude  of  the  backlash  angle.  This  fact  is  significant  since  it 
permitted  the  use  of  a  constant  backlash  angle  in  the  solution  of  the 
problems,  yet,  in  effect,  produced  results  for  any  amount  of  backlash. 
Secondly,  the  amplitude  of  maximum  overshoot  at  zero  backlash  angle 
agrees  with  that  obtained  by  the  analytical  solution  of  the  second 
order  linear  system.  This  agreement  is  taken  as  a  verification  of  the 
accuracy  of  the  computer  solutions. 
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Chapter  V 
DISCUSSION  OF  RESPONSE  CHARACTERISTICS 

A.  Maximum  Overshoot  for  the  Plastic  Impact  Case. 

After  numerous  data  had  been  accumulated,  several  schemes  of  pre- 
senting the  information  were  attempted  to  find  an  optimum  method  for 
use„  One  such  attempt  is  shovm  in  figure  8  which  is  a  plot  of  maxinium 
overshoot  versus  the  ratio  of  load  friction  to  the  total  friction  at  a 
constant  damping  coefficient.  This  method  was  found  to  have  serious 
disadvantages.  Because  of  appreciable  crossing  of  the  lines  of  con- 
stant inertia  ratio^  it  was  found  difficult  to  interpolate  throughout 
the  desired  range  of  inertia  ratios.  Also  there  was  generally  a  small 
variation  of  maximum  overshoot  with  friction  ratios. 

Another  method  which  was  attempted  was  a  plot  of  maximum  over- 
shoot versus  damping  coefficient  at  constant  friction  and  inertia 
ratios.  These  curves  showed  the  effect  of  various  inertia  and  fric- 
tion ratios  in  a  system  with  backlash  on  the  linear  system  curve  of 
overshoot  versus  damping  coefficient  (figure  9),  but  was  foxind  diffi- 
cult to  use.  The  nianber  of  curves  needed  to  present  the  data  would 
have  been  excessive.  Also  interpolation  throughout  the  desired  range 
of  inertia  ratios  would  have  been  extremely  difficult. 

The  curves  of  figures  10  through  15,  which  are  plots  of  maximum 
overshoot  versus  J,/J^,  were  found  to  be  the  most  desirable  method  of 
presenting  the  data.  None  of  the  disadvantages  apparent  in  the  other 
methods  is  present  here.  From  these  six  charts  an  accurate  prediction 
of  the  maximum  overshoot  can  be  made  for  any  arbitrary  second  order 
servo  system  \-jith  backlash  using  the  following  equations 
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The  use  of  these  charts  in  predicting  system  response  is  explained  and 
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-   M 


illustrated  in  Appendix  A< 
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As  was  pointed  out  earlier^  the  solutions  were  made  and  th|'  results 
were  plotted  using  a  backlash  angle  of  0.3  radian.  Thus  the  values  of 
amplitude  of  maximum  overshoot  on  the  curves  are  significant  only  in  that 
thej'  represent  the  slope  of  the  overshoot  versus  backlash  straight  lines. 
This  slope  must  be  calculated  in  order  to  scale  the  value  of  overshoot 
read  from  the  curve  for  the  magnitude  of  backlash  actually  expected. 
Therefore^  in  order  to  refine  the  prediction  scheme^  other  design  curves^, 
figures  16  through  21,  were  calculated  using  the  curves  of  figures  10 

dl'i 

through  15  just  discussed.  Along  the  curves  the  slope   pt 

calculated  by^ 


was 


vX 


d  A 


f'f. 


pt 


-   M 


p^Ia-o 


4- .3 


(V^2) 


These  values  of  dM  ^/dijk   were  then  plotted  versus  inertia  ratios  at 

pt' 

constant  friction  ratios  as  before.  These  curves  give  a  more  straight- 
forvrard  pre'^iction  method  with  the  following  equations 


pt 


pt  pt' 

^=0  1/=;^ 


X 


A, 


(V-3) 


The  use  of  these  design  charts  is  explained  -and  illustrated  in  Appendix  A 
along  vd-th  the  former  nethod. 

Figure  23  shoifs  a  graphical  presentation  of  the  transient  response 
of  the  system  containing  various  parfimeters.  A  comparison  of  curves  (3) 
and  (4)  readily  shows  the  effect  of  backlash  on  the  maximum  overshoot 
with  ether  parameters  remaining  constant.  A  comparison  of  curves  (1) 
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and  (4)  showss  the  effect  of  varying  the  damping  ratio  of  the  system  \i±ih 
other  paiameters  remaining  constant. 
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Bo  Maxiaom  Overshoot  for  the  Elastic  Impact  Case. 

One  of  the  assumptions  made  in  the  beginning  of  this  study  was 
that  the  contact  of  the  gear  teeth  is  plastic  with  no  bouncing  of  the 
gears o  Since  this  study  was  begun ^  an  investigation  has  been  completed 
at  the  Uo  So  Postgraduate  School  which  considers  the  backlash  case  with 
varyiiTig  amounts  of  resilience  in  the  gear  teeth.  This  investigation^, 
carried  on  by  T„  Wo  Luckettj,  and  No  Oo  Anderson  has  produced  some  resulti 
bearing  on  max.lroujrt  overshoot  which  will  be  included  here.  The  informa= 
tion^  while  not.  included  as  part  of  their  study,  was  available  from  the 
data  taken  from  the  computer  solution  of  their  problem. 

The  analysis  of  this  data  revealed  a  surprising  effect  of  gear 

boiinc«  on  tne  maximi'm  overshoot  of  the  system.  The  amount  of  over= 

shoot  with  plastic  impact  is  always  greater  than  that  with  elastic 

impact o  The  difference  between  this  maximum  overshoot  decreases  as 

F,/F.  ^  ,  i^icreaseso  However^  under  no  circumstances  was  the  effect 
L  total  ^ 

on  maximuir  overshoot  greater  than  .01  radians.  The  effect  on  the 
time  of  rria:<im'am  overshoot  appears  negligible.  In  a  few  cases  the  time 
was  different  for  the  plastic  and  the  elastic  cases  but  the  difference 
was  very  small  and  no  pattern  of  difference  could  be  ascertained. 

These  results  point  out  the  validity  of  the  original  assumption 
made  for  this  study o  In  a  practical  system  there  will  exist,  of 
course^,  some  amount  of  gear  boimce  v/hen  the  gears  contact  each  other. 
The  effect  en  maximim  overshoot  ^  however,  is  so  small  that  the  curves 
obtained  in  this  study  remain  reasonably  accurate. 
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Co  SettliniR  Time  and  Time  of  Maximum  Overshoot , 
lo  Settljjig  TiSl? 

In  deterrdning  the  effect  of  backlash  on  the  settling  time  of  a 
second  order  servo  systemi,  it  was  necessary  to  first  find  a  definition 
for  settling  time   which  would  be  suitable  for  all  of  the  cases  which 
were  investigated.  Normally ^  settling  time  is  defined  as  the  time  for 
the  output  to  reach  9B^  of   its  final  value.  This  definition  has  real 
meaning^,  however^  only  for  linear  stable  systems  in  which  there  are  no 
sustained  oscillationso  For  systems  which  do  exibit  small  sustained 
oscillatioiicy,  the  settling  time  is  sometimes  defined  as  the  time 
required  for  the  amplitude  of  oscillations  to  decrease  to  within  2% 
of  the  final  value o  For  very  small  amplitudes  of  oscillations  this 
definition  gives  an  unrealistic  value  for  settling  time.  For  the 
purpo&h   of  this  study  settling  time  is  considered  as  the  time  required 
for  the  system  output  to  settle  to  ^S%   of  the  difference  between  the 
magnitude  of  the  input  and  the  amplitude  of  the  sustained  oscillations 
as  shown  in  the  sketch  below. 


6u 


TiKE.  ('SECOUDS) 
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For  all  cases  for  which  the  damping  coefficient  was  less  than  0.6 
there  were  no  conclusions  as  to  effects  of  backlash  on  the  system 
settling  time  since  the  magnitudes  of  the  sustained  oscillations  were 
always  quite  large  for  the  amount  of  backlash  used  in  the  solutions  and 
the  settling  time  was  difficult  to  obtain  accurately.  For  those  cases 
for  which  the  damping  coefficient  was  greater  than  0.6,  the  sustained 
oscillations  were  small  if  they  existed  and  an  analysis  revealed  the 
following  general  conclusions s 

lo  As  the  size  of  the  backlash  angle  increases,  the  settling  time 
will  increase  provided  the  other  parameters  are  kept  constant o  For  the 
cases  investigated^  this  increase  in  settling  time  was  nearly  proportion- 
al  to  the  backlash  angle, 

2,  As  the  ratio  ^-r/^-^   increases,  with  other  parameters  constant^ 
the  settling  time  increases, 

3.  At  the  ratio  F./F.  x  nis  increased,  with  other  parameters 
constant^  the  settling  time  decreases. 

The  data  included  in  Tables  I-III  show  these  general  conclusions 
and  also  some  q\iantitative  effects  of  backlash  on  settling  time. 
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TABLE  I 
SETTLING  TBIE  DATA 


/ 

total 

Radians 

A 

Seconds 

t 

s 

06 

1/9 

.2 

.3 

28 

06 

1/9 

o4 

.3 

21 

06 

1/9 

.6 

.3 

16 

»6 

1/9 

.8 

.3 

13 

„6 

1/9 

1.0 

.3 

7 

,6 

1/4 

,2 

.3 

28 

06 

1/4 

o4 

.3 

24 

.6 

1/4 

.6 

.3 

16 

,6 

1/4 

.Q 

0 

U 

06 

1/4 

1.0 

.3 

9 

e6 

1 

.2 

.3 

29 

.6 

1 

.4 

.3 

28 

.6 

1 

.6 

.3 

23 

.6 

1 

.a 

.3 

17 

.6 

1 

l.C 

.3 

10 

,6 

4 

.2 

.3 

29 

.6 

4 

o4 

.3 

28 

06 

4 

.6 

.3 

26 

.6 

4 

.a 

0 
•-> 

23 

.6 

4 

loO 

.3 

10 

.6 

9 

,2 

.3 

29 

.6 

9 

.4 

.3 

29 

.6 

9 

06 

.3 

26 

.6 

9 

.a 

.3 

23 

,6 

9 

1.0 

.3 

10 
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TABLE  II 
SETTLING  TIME  DATA 


/ 

'^1 

/ 
^ total 

Radians 

Seconds 

t 

s 

.6 

1/9 

.8 

.03 

5 

.6 

1/9 

.3 

.10 

7 

06 

1/9 

c8 

.15 

8 

.6 

1/9 

08 

.30 

13 

,6 

1 

.8 

,03 

11 

06 

J. 

08 

.10 

U 

.6 

„8 

as 

16 

06 

1 

,8 

.30 

17 

e6 

,8 

.03 

12 

c6 

k 

.8 

.10 

lA 

.6 

V 

.8 

0I5 

17 

c6 

4 

,8 

.30 

23 

9 

.8 

oC3 

13 

.6 

9 

.8 

.10 

16 

.6 

9 

.8 

.15 

18 

.6 

9 

.8 

.30 

23 
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TABLE  III 
SETTLING  TBE  DATA 


y 

V 
'  F 
total 

Radians 

A 

Seconds 

4- 

s 

.8 

1/9 

.2 

.3 

8 

.8 

1/9 

.4 

.3 

7 

.3 

1/9 

.6 

.3 

6 

.8 

1/9 

.8 

.3 

5 

.8 

1/9 

1.0 

.3 

5 

,8 

1/4 

.2 

.3 

20 

.8 

1/4 

.4 

.3 

9 

o8 

1/4 

.6 

.3 

7 

,8 

1/4 

.8 

.3 

6 

,8 

1/4 

1.0 

.3 

5 

.8 

1 

.2 

.3 

43 

o8 

1 

.4 

.3 

31 

.8 

1 

.6 

.3 

10 

.8 

1 

.8 

.3 

8 

o8 

1 

1.0 

.3 

6 

,8 

4 

.2 

.3 

48 

,8 

4 

.4 

.3 

36 

,8 

4 

.6 

.3 

12 

.8 

4 

.8 

.3 

9 

,8 

4 

1.0 

.3 

6 

c8 

9 

.2 

.3 

52 

.8 

9 

.4 

.3 

43 

.8 

9 

.6 

.3 

13 

o8 

9 

.8 

.3 

9 

.8 

9 

1.0 

i. 

.3 

6 
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^»  Time  of  Maximum  Overshoot 

For  each  of  the  cases  investigated  the  time  of  maximum  overshoot 
(peak  time)  was  recorded  from  the  printout  of  the  solution.  It  was 
found  that  the  presence  of  backlash  causes  a  slight  increase  in  the 
peak  time  over  the  linear  system  for  those  cases  where   JT  ^.5^  and 
has  negligible  effect  in  those  where   J*^  o5  »  When  Fj/Fx  ^  n  is 
increased  from  zero  to  0.5,  other  parameters  constant  and  J  x^  .5, 
peak  time  increases  slightly  and  remains  essentially  constant  for 
F-/F,  .  1  ^  .5  .  As  Jr/Jiyr  IS  Increased,  other  parameters  constant 
and   J  X*  .5>  peak  time  is  decreased o 


U2 


D,   Other  System  Characteristics . 

A  careful  study  of  the  results  as  plotted  in  figures  10-15  reveals 
some  interesting  phenomena.  At  all  friction  ratios  it  is  observed  that 
as  the  ratio  of  load  inertia  to  motor  inertia  is  made  smaller  the  ampli- 
tude of  tlie  overshoot  is  less.  There  is  a  point  where  the  overshoot  is 
identical  to  that  of  a  linear  second  order  system^  in  other  words,  where 
backlash  has  no  effect  on  the  overshoot.  Note  that  this  point  occurs 
at  different  inertia  ratios  for  each  value  of  damping  ratio.  Also  this 
point  varies  depending  on  the  value  of  friction  ratio.  If  the  inertia 
ratio  is  decreased  beyond  this  point  then  the  amplitude  of  overshoot 
remains  constant  at  the  value  of  the  linear  case.  Figure  22  is  a  plot 
of  the  boundaries  of  the  combinations  of  variables  which  will  cause  an 
effect  on  the  araplitxide  of  overshoot  if  tlie  system  has  any  appreciable 
amount  of  backlash.  The  coordinates  are  friction  ratio  and  inertia 
ratio  with  the  boundaries  labeled  with  values  of  damping  ratio.  Areas 
to  the  right  of  a  boundary  are  combinations  of  the  variables  which 
affect  overshoot.  These  curves  were  plotted  using  inforniation  deducted 
from  figures  10-15 «  Three  points  were  selected  as  test  points  to  prove 
the  validity  of  the  boundaries  as  plotted.  These  points  are  shown  on 
the  figure  as  points  "A",  "B",  and  "C".  The  system  parameters  used  at 
the  test  points  were? 

Point  A  Point  B  Point  C 

/     =0.8       y     =  0.6       f  ^   0.3 

V^total     =°-^  V^total     =^'^  V^total     ==  ^'^ 
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The  system  was  solved  or.   the  digital  computer  to  obtain  the  raaxinujn 
overshoot  for  these  three  sets  o£   parameters.  According  to  the  ciirves 
point  A  falls  to  the  left  of  the  boundary  of  the  J  =  .3  area. 
Therefore,  the  overshoot  should  be  the  same  as  that  of  the  linear  case. 
The  computer  solution  gave  an  overshoot  of  1.019,  which  is,  in  fact, 
identical  to  that  of  the  linear  case  as  seen  in  figure  9«  Similarly, 
point  B  falls  to  the  right  of  the  J  =  0.6  boundary  and  in  the  area 
where  the  amplitude  of  overshoot  should  not  be  th.e  same  as  for  the 
linear  case.  The  computer  solution  gave  a  maximum  overshoot  of  1.110 
which  is  slightly  greater  than  the  linear  case.  Point  G  falls  to  the 
right  of  the  J  =  0,3  boundary  and  in  the  area  where  the  amplitude  of 
overshoot  should  be  different  from  that  of  the  linear  system.  The 
computer  solution  gave  a  maximum  overshoot  of  l.Al  while  that  of  the 
linear  system  is  1,372.  From  these  boundary  curves  one  cannot  deter- 
mine the  exact  effect  of  backlash  on  the  maxifnum  overshoot  ^f  a  system, 
but  can  deterriiine  in  which  systems  backlash  will  cause  an  appreciable 
effect. 

Another  phenomena  which  is  revealed  ir.  the  results  as  plotted  in 
figui-es  10-15  is  thi,t  for  all  cases  the  greatest  maximum  overshoot 
occurs  in  the  range  0,1  \   «^t /^^^^  \  1^»0  ^^^  ^^^-   rarigo  of  Jj/'^y,  test- 
ed.    However,  it  appears  that  since  the  curves  are  trending  upwards  at 
the  higli  end  of  the  Jt/«J»t  axis  that  at  very  large  inertia  ratios  the 
overshoot  would  be  greater.  It  is  felt  that  these  large  inertia  ratios 
are  outside  the  practical  range.  An  attempt  v;as  made  to  explain  these 
phenomena  but  without  success.  It  is  often  very  difficult,  if  not 
impossible^ to  find  a  physical  explanation  for  certain  behavior  caused 
by  non-linearities.  They  can  merely  be  pointed  out  from  analyses  of 
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experimental  information  and  accounted  for  accordingly. 

It  is  interesting  to  note  that  when  the  ratio  of  load  friction  to 
motor  friction  is  below  0.4,  as  in  figures  10  and  11,  there  is  a  possi- 
bility of  an  overshoot  for  the  critically  damped  case,  i.e.  J  =  1.0, 
depending  on  the  inertia  ratio.  At  first  glance  one  would  reason  that 
if  the  system  is  critically  damped  there  should  be  no  overshoot.  How- 
ever, it  can  be  seen  that  if  the  friction  ratio  is  small  and  the  inertia 
ratio  is  large  then  the  ratio  of  F, /j, ,  which  is  the  slope  of  the  drift- 
ing load  trajectory,  is  made  small.  Therefore,  when  the  load  separates 
from  the  motor,  its  velocity  remains  fairly  constant  at  a  relatively 
high  value  and  its  momentum  carries  it  over  the  input  to  cause  an 
overshoot. 
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Chapter  VI 
CONCLUSIONS 

The  following  is  a  siommary  of  facts  which  can  be  concluded  from  the 
results  of  this  investigation: 

1.  The  amplitude  of  mai^cimum  overshoot  in  a  second  order  servo  sys- 
tem with  backlash  bears  a  linear  relationship  with  the  magnitude  of  the 
backlash  angle.  The  slope  of  the  straight  line  relating  overshoot  to 
the  backlash  angle  varies  depending  on  the  division  of  inertia  betv;een 
motor  and  load,  the  division  of  friction  between  motor  and  load,  and  the 
system  damping  coefficient.  The  design  charts  included  in  this  thesis, 
drawn  from  the  results  of  many  solutions  of  the  equations  of  the  motion 
of  the  system  and  based  on  a  backlash  angle  of  0.3  radian^  can  be  used 
to  predict  the  amount  of  overshoot  of  a  particular  servo  system. 

2.  The  fact  that  the  gear  teeth  have  a  certain  amount  of  resilience 
and  the  contact  is  not  plastic,  has  negligible  effect  on  the  results  ob- 
tained in  this  study, 

3.  The  primary  effect  of  backlash  on  the  settling  time  of  a  second 
order  servo  is  to  increase  it,  the  increase  being  proportional  to  the 
magnitude  of  the  backlash  angle.  Due  to  the  presence  of  limit  cycles 

of  sizable  magnitudes  for  systems  with  -^    <.6,  settling  time  for  those 
systems  was  difficult  to  define  and 'no  conclusions  v;ere  reached.  For 
those  systems  with  -j  y   .6?  (a)  Increasing  ^-J'^^^^^i   decreases  settling 
time,  (b)  Increasing  Jr/^vr  increases  settling  time,  (c)  Increasing  the 
damping  coefficient  increases  settling  time,  and  (d)  Increasing  back'^ 
lash  increases  settling  time, 

4.  The  presence  of  backlash  in  a  system  causes  a  slight  increase 

in  the  tine  of  maximum  overshoot  (peak  time)  for  those  systems  where^>,5i 

IS 


and  has  negligible  effect  in  those  where  T'<'o5o  When  F,/F^  .  ,  is 

-'  L     total 

increased  from  zero  to  0.5,  other  parameters  constant  and  jf]^  .5j. 
peak  time  increases  slightly  and  remains  essentially  constant  for 
F^/F.  + , -]  >  .5.  As  ^j/J^   is  increased,  other  parameters  constant  and 
y  ^   *^i   peal<:  time  is  decreased, 

5.  With  certain  combinations  of  parameters  present  in  a  system, 
backlash  effects  no  change  in  the  amplitude  of  overshoot  over  that  of 
the  linear  system.  Regions  of  these  combinations  are  jjresented  in  chart 
form  in  figure  22. 

6.  For  all  friction  ratios  and  damping  coefficients  the  greatest 
maximum  overshoot,  when  backlash  is  present,  occ\irs  in  the  range 

0.1  <  Jj./J„<  10. C  for  the  range  of  «Jt/J^j  investigated. 

7.  In  systems  where  F,/F.  ^  -,  <^  .4,  when  backlash  is  present, 
there  will  be  an  overshoot  even  for  critically  damped  systems  if  ^j/\i! 

is  high  enough. 
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APPENDIX  A 
USE  OF  DESIGN  CHARTS 

In  order  to  demonstrate  the  procedure  for  using  the  design  charts 
shovm  in  figures  10-21,  an  example  problem  is  solved  obtaining  the 

maximum  overshoot,  M  . ,  for  a  second  order  servo  with  backlash.  From 

'     pt' 

Ref .   2,   the  differential  equation  for  the  combined  motion  of  an  armature 
controlled  motor  is: 

Q    +     £kjLlL.Q    ^        ^       Q  ^      ^       Q   (A-l) 


Given  the  following  specifications  of  a  D.  C.  shunt  motor  coupled 
through  a  gear  box  to  a  load: 
Horsepower  =  .008 
Rated  RPM   =  4000 
R  =  282  ohms 

K         =  .233  volts/radian/second 

K.         =  .1464  lb . =f eet/ampere 
J         =  8.2  X  10"^  slug-feet^ 


m 

F 
m 


=  3.34  X  10~  lb. -feet/radian/sec. 


K  =10  volts/radian 

e 


^ 


=  1 


A       =  .05 


J         =  16.4  X  10"^  slug-feet^ 


^L 


=  10.4  X  10   lb. -feet /radian/sec. 
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The  required  values  shown  in  Equation  A-1  are  obtained  as  follows; 


K 
a 


~6 


—I 


,-     ^/qo 


-6 


/3.7^    ^/0~^      ^^--^^ 


^^^/sec 


^AD/S£c 


"^L/^otal    /J.J-  X  /^  -^     •  ^^7 


Two  methods  are  shown  for  obtaining  the  maxiraura  overshoot  with  a 
backlash  angle  of  .05  radians  in  the  gear  trains. 

The  first  method  is  accomplished  using  the  curves  of  figures  10-15. 
The  procedure  is  simple,  but  requires  more  steps  than  does  the  second. 
Also  one  must  remember  that  these  curves  are  drawn  for  a  backlash  of 
0.3  radian  and  must  be  scaled  for  the  expected  amount  of  backlash.  The 
second  method  makes  use  of  the  design  charts  of  figures  16-21. 
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Method  1 


Entering  figure  10,  (Ft/F+  +  1  ~  '^'■^^»   ^^   '^L^'^M  ~  '^"^^   ^^^  ^   danip= 
ing  ratio, j      =  0,19,  the  amount  of  overshoot  for  a  backlash  angle^  ^  - 


0.3  rad.  is  M  .  -  1.658. 
pt 


Using  Equation  V-1 


M 


•pt 


f=..l9 


=  M        +  .05/.3 


M 


pt 


M 


pt 


_) 


From  the  second  order  linear  characteristics  curve  of  figure  9 


hence. 


V  ^  ^-^"^^ 


Mp^  =  1.545  +  .05/.3 
=  1.545  +  .0189 


1.658 


1.545 


Method  2 


Using  Equation  V-3 


M 


•pt 


M  ^     +   dM  ./dA 


Z^ 


f^JB 


from  figure  16,  (FLAtotal  "  ^'"^^^  ^^  *^l/''^M  "  2.0,  and  a  damping  ratio, 


=  0.19,  we  obtain    dM  ^/dA  =  0,357 


From  figure  9 


M 


■pt 


=  1.545 


hence. 


M  ^  =  1.545  +  .357  X  .05 
=  1.545  +  .018 
=  1»563 
This  compares  closely  to  the  answer  obtained  using  Method  1. 
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The  procedure  for  using  the  charts  remains  unchanged  for  servo 
systems  which  utilize  other  types  of  motors.  The  only  difference  which 
arises  is  in  the  calculation  of  the  chart  pararaeters  from  the  knoim 
system  constants.  For  a  field  control  motor  if  the  field  time  constant 
is  negligible,  the  equation  is  similar  to  the  armature  control  motor 
equation. 


9  f  -r^  9  f  ^^  0  -  y^  4     (A-2) 


In  like  manner  the  equation  for  the  system  with  a  two  phase  motor  is 


where 


\\       =  -  ^ —    (from  characteristic  curves  of  two  phase  motor) 


K 


=  Ol  (from  characteristic  curves  of  two  phase  motor) 


m     -N 


The  calculation  of  the  chart  parameters  is  a  simple  matter  once  the 
defining  equations  are  known. 
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APPEIJDIX  B 


COMPUTER  PROGRAM 


REM  BACKLASH  PR0BLEI«4  PLASTIC 
CONTACT  RUN  3 


ORG 

20000 

20000 

75  4  20022 

START 

SLJ 

4 

njpuT  1 

REMARKS  IN  THIS 

00  0  00000 

0 

0 

0 

COLUMN 

20001 

75  4  20063 

COMBINED 

SLJ 

4 

RESET 

00  0  00000 

0 

0 

0 

20002 

75  4  60200 

SLJ 

4 

RUNGE 

SET  UP  GILL 

50  0  00000 

ENI 

0 

0 

(Procedure  to  solve 

20003 

00  0  20023 

0 

0 

TABLE  1 

differential 

00  0  20034 

0 

0 

DERC 

equations) 

200C4 

75  4  60201 

POINTC 

SLJ 

4 

RUNGE  +  1 

00  0  00000 

0 

0 

0 

20005 

75  4  60200 

LOAD  1 

SLJ 

4 

RUNGE 

SET  UP  GILL 

50  0  00000 

ENI 

0 

0 

20006 

00  0  20040 

0 

0 

TABLE  L 

00  0  20051 

0 

0 

DERL 

20007 

75  4  60201 

POINTL 

SLJ 

4 

RUNGE  +  1 

00  0  OOfjOO 

0 

0 

0 

20010 

75  4  20130 

SLJ 

4 

COWJT 

50  0  00000 

ENI 

0 

0 

20011 

75  0  20061 

SLJ 

0 

TEST  1 

50  0  00000 

ENI 

0 

0 

20012 

75  4  60200 

MOTOR 

SLJ 

4 

RUNGE 

SET  UP  GTTT. 

50  0  00000 

ETJI 

0 

0 

20013 

00  0  20023 

0 

0 

TABLE  1 

00  0  20054 

0 

0 

DERM 

20014 

75  4  60201 

POINTM 

SLJ 

4 

RUNGE  +  1 

00  0  00000 

0 

0 

0 

20015 

75  4  60200 

L0AD2 

SLJ 

4 

RUI\IGE 

SET  UP  GILL 

50  0  00000 

sr:i 

0 

0 

20016 

00  0  20040 

c 

0 

TABLE  L 

00  0  20051 

0 

0 

DERL 

20017 

75  4  60201 

P0INT2L 

su 

4 

RUNGE  +  1 

OOC  0  00000 

0 

0 

0 

200-^.0 

75  4  20130 

SLJ 

4 

COUNT 

50  0  00000 

ENI 

0 

0 

20021 

75  0  20073 

SLJ 

0 

TEST  2 

00  0  00000 

0 

0 

0 

20022 

75  0  00000 

INPUTl 

SLJ 

0 

0 

Provision  for  chang- 

75 0  20022 

SLJ 

0 

INPUTl 

ing  Program  if  desire^ 

20023 

00  0  00000 

TABLEl 

DEC 

2 

•Jable  for  gill  combii 

00  0  00002 

20024 

17  7  15075 
34  1  21727 

DEC 

.01 

20025 

00  0  00000 
00  0  00000 

T 

DEC 

0 

56 


20026 

00 

00 

0 
0 

00000 
00000 

UDOT 

DEC 

0 

20027 

00 
00 

0 
0 

00000 
00000 

U 

DEC 

0 

20030 

00 
00 

0 
0 

00000 
00000 

QU 

DEC 

0 

20031 

00 
00 

0 
0 

00000 
00000 

THETADOT 

DEC 

0 

20032 

00 

0  "> 

0 
0 

00000 
00000 

THETA 

DEC 

0 

20033 

00 

00 

0 

0 

00000 
00000 

QTHETA 

DEC 

0 

20034 

13 

0 

20156 

DERC 

LAC 

0 

3          Derivative  for  com- 

32 

0 

20027 

F^aJ 

0 

y          bined 

20035 

31  0 

20032 

FSB 

0 

THETA 

30 

0 

20164 

FAD 

0 

ONE 

20036 

20 

0 

20026 

STA 

0 

UDOT 

12 

0 

20027 

LDA 

0 

U 

20037 

20 

0 

20031 

STA 

0 

THETADOT 

75 

0 

60202 

SLJ 

0 

RUNGE  +  2 

20040 

00 
00 

0 
0 

00000 
00002 

TABLET, 

DEC 

2          Table  for  gill  loac 

20041 

17  7  15075 
34  1  21727 

DSC 

.01 

20042 

00 
00 

0 
0 

00000 
00000 

TL 

DEC 

0 

20043 

00 
00 

0 

0 

00000 
00000 

vrooT 

DEC 

0 

20044 

00 

00 

0 
0 

00000 
00000 

w 

DEC 

0 

20045 

00 
00 

0 

00000 
00000 

aw 

DEC 

0 

20046 

00 
00 

0 
0 

00000 
00000 

TPffiTAT,n 

DEC 

0 

20047 

00 
00 

0 
0 

00000 
00000 

TIIETAL 

DEC 

0 

20050 

00 
00 

0 
0 

00000 
00000 

QTHETAL 

DEC 

0 

20051 

13 

0 

20157 

DEIiL 

LAG 

0 

D          Derivative  for  loac 

32 

0 

20044 

n'!U 

0 

W 

20052 

20 

0 

20043 

STA 

0 

i-rooT 

12  0  20044 

LDA 

0 

W 

20053 

20  0  20046 

STA 

0 

THETA  T.n 

75 

0 

60202 

SLJ 

0 

RUNGE  +  2 

20054 

13 

0  20161 

DSRI'-I 

LAC 

0 

F          Derivative  for  mctc 

32  0  20027 

F14U 

0 

U 

20055 

31 

0 

20047 

FSB 

0 

THETAL 

30  0  20164 

FAD 

0 

ONE 

20056 

33 

0 

20160 

FDV 

0 

E 

20 

0 

20026 

STA 

0 

UDOT 

20057 

12 

0 

20027 

LDA 

0 

U 

20  0  20031 

STA 

0 

THETADOT 

57 


20060 

75 

0 

60202 

SLJ 

0 

RUNGS  +  2 

50  0 

00000 

ENI 

0 

0 

20061 

12 

0 

20046 

TESTl 

LDA 

0 

THETHAID 

Test  to  determine 

31  0 

20031 

FSB 

0 

TKETaDOT 

Start  cf  backlash 

20062 

22 

2 

20012 

JUMFl 

AJP 

2 

MOTOR 

75 

0 

20001 

SLJ 

0 

COlffilNED 

20063 

75 

0 

00000 

liESET 

SLJ 

0 

0 

12 

0 

;-:0025 

LDA 

0 

T 

Changes  load  initial 

20064 

20 

0 

20042 

STA 

0 

TL 

conditions  back  to 

12 

0 

20026 

LDA 

0 

UDOT 

motor  initial  con- 

20065 

20 

0 

20043 

STA 

0 

IVDOT 

ditions  if  separa- 

12 

0 

20027 

LDA 

0 

U 

tion  has  not  occurr- 

20066 

20 

0 

20044 

STA 

0 

\-i 

ed 

12 

0 

20031 

LDA 

0 

tpiet;j)ot 

20067 

20 

0 

20C46 

STA 

0 

THETALD 

12 

0 

20032 

LDA 

0 

THETA 

20070 

20 

0 

20047 

STA 

0 

THETAL 

10 

0 

00000 

SljUti 

0 

0 

20071 

20 

0 

20045 

STA 

0 

QW 

20 

0 

20050 

STA 

0 

QTHETAL 

20072 

75  0  20063 

SLJ 

0 

RESET 

50 

0 

00000 

EfJI 

0 

0 

20073 

12 

0 

20047 

TEST2 

LDA 

0 

THET.y: 

Test  to  determine  if 

31 

0 

20032 

FSB 

0 

THETA 

backlash  has  be  ^n 

20074 

31  0 

20166 

BKTEST2 

FSB 

0 

BACKLASH 

taken  up 

50 

0 

00000 

EIII 

0 

0 

20075 

22 

2 

20137 

jni-2 

AJF 

o 

COOTACT 

75 

0 

20012 

SLJ 

0 

MOTOR 

20076 

36  0  20117 

SVvTTOH 

SSK 

0 

ALTERNAT 

Plan  for  changing 

75  0  20100 

SLJ 

0 

LEFT 

sign  of  backlash 

20077 

75 

0 

20107 

SLJ 

0 

RIGHT 

each  cycle 

50 

0 

00000 

ENI 

0 

0 

20100 

12 

0 

20062 

LEFT 

LDA 

0 

JUMP  1 

14  0 

20116 

ADD 

0 

FACTOR 

20101 

20 

0 

20062 

STA 

0 

JUMP  1 

12 

0 

20075 

LDA 

0 

JUI-IP  2 

20102 

14  0 

20116 

;uDD 

0 

FACTOR 

20 

0 

20075 

STA 

0 

JUMP  2 

20103 

12 

0 

20074 

LDA 

0 

3KTE3T2 

15 

0 

20172 

SUB 

0 

FACT0R2 

20104 

20  0  20074 

STA 

0 

BKTEST2 

12  0  20117 

LDA 

0 

ALTERI-.^AT 

20105 

05 

0 

00001 

AJ3 

0 

1 

20 

0 

20117 

STA 

0 

ALTERIjAT 

20106 

75 

0 

20001 

SLJ 

0 

COMBINED 

50 

0 

00000 

EMI 

c 

V^' 

20107 

12 

0 

20062 

RIGHT 

LDA 

0 

JUMPl 

15 

0 

20116 

SUB 

0 

FACTOR 

20110 

20 

0 

20062 

STA 

0 

JUl^IPl 

12 

0 

20075 

LDA 

0 

jmT2 

20111 

15  0 

20116 

SUB 

0 

FACTOR 

J 1  , 

0 

20075 

Jin 

0 

JU1-IP2 

20112 

12 

0 

20074 

LDA 

0 

BKTSST2 

58 


14  0  20172 

/iDD 

0 

FACTCR2 

20113 

20  0  20074 

STA 

0 

BKTEST2 

12  0  20117 

LJA 

0 

fiJJIKi'AAT 

20114 

05  0  00001 

AL3 

0 

1 

20  0  20117 

STA 

0 

ALTERNAT 

20115 

75  0  20001 

SLJ 

0 

COUBJNED 

50  0  00000 

ENI 

0 

0 

20116 

00  1  00000 

FACTOR 

OCT 

0010000000000000 

00  0  oocoo 

20117 

25  2  52525 

ALTEiiNAT 

OCT 

2525252525252525 

25  2  52525 

20120 

75  0  00000 

FdJJiiT 

SLJ 

0 

0 

Print  routine 

12  0  20025 

IDA 

0 

Tl 
X 

20121 

20  0  20252 

STA 

0 

BUF 

12  0  20027 

LDA 

0 

U 

Motor  speed 

20122 

20  0  20253 

STA 

0 

BUF  +  1 

12  0  20032 
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